Labyrinthine clustering in a spatial rock-paper-scissors ecosystem 



Jeppe Juul, Kim Sneppen, and Joachim Mathiesen 
University of Copenhagen, Niels Bohr Institute, Blegdamsvej 17, DK-2100 Copenhagen, Denmark 

(Dated: October 29, 2012) 

The spatial rock-paper-scissors ecosystem, where three species interact cyclically, is a model ex- 
ample of how spatial structure can maintain biodiversity. We here consider such a system for a 
broad range of interaction rates. When one species grows very slowly, this species and its prey 
dominate the system by self-organizing into a labyrinthine configuration in which the third species 
propagates. The cluster size distributions of the two dominating species have heavy tails and the 
configuration is stabilized through a complex, spatial feedback loop. We introduce a new statis- 
tical measure that quantifies the amount of clustering in the spatial system by comparison with 
its mean field approximation. Hereby, we are able to quantitatively explain how the labyrinthine 
configuration slows down the dynamics and stabilizes the system. 

PACS numbers: 87.23.Kg, 87.23.Cc, 87.18.Hf 



Introduction — Spatial migration of species is crucial 
for the viability of many ecological systems. As a striking 
example, crickets are known to locally deplete their nutri- 
tional resources to an extend where mass-migration is the 
only alternative to cannibalism pii i2j . Once the crickets 
have left an area, they can not return until the natu- 
ral resources have been reestablished. Likewise, deadly 
viruses and bacteria depend on constantly infecting new 
hosts to survive [3-5 . 

The rock-paper scissors game has emerged as a 
paradigm to describe the impact of spatial structure on 
biodiversity [6-11 . In this system, three species interact 
cyclically such that species 1 can invade species 2, which 
can invade species 3, which, in turn, can invade species 
1 (see Fig. [l^). Such intransitive interaction pattern is 
very similar to the important genetic regulatory network 
the repressilator [12", T3^ and has been identified in many 
ecological system, among others in marine benthic sys- 
tems [m [15!, plant systems [TBHIH], terrestial systems 
[I9l|20], and microbial systems [2TH25] . In such systems, 
all species constantly need to migrate spatially to sur- 
vive. Investigating three strands of E. coli bacteria with 
cyclic interactions, it has been shown that biodiversity 
can not be preserved unless spacial structure is imposed 
by arranging the bacteria on a petri dish [71 [HI [26] . These 
results have been reproduced in Monte Carlo simulations 
[6^ 27-29 , but even though many different analytical ap- 
proaches have been applied, exactly how spatial structure 
stabilizes the system is still an open problem [3QH33] . 

Model — We study the rock-paper-scissors game on 
a square lattice of L x L nodes and periodic boundary 
conditions. Each node is occupied by one of the three 
species 1, 2, or 3 growing at rates Vi^ and vs^ respec- 
tively. In each update a random node i and a random of 
its neighbors j are selected. If i can invade j according 
to the cyclic interacting pattern illustrated in Fig. [T^, it 
will do so with a probability equal to Vi. 

Results — When the three species are initiated from a 
random configuration and with equal growth rates, they 



quickly organize into a steady state where all species are 
equally abundant and form small clusters (see Fig. ^p)- 
If the growth rate of species 3 is increased compared to 
species 1 and 2, species 2 becomes more abundant on the 
lattice and all three species form larger clusters (see Fig. 
[l]^). This paradoxical behavior, that the biomass of one 
species increases proportional to the growth rate of its 
prey, is characteristic for the rock-paper-scissors system 

El Hz]. 

Similarly, if the growth rate of species 1 is decreased, 
species 3 slowly becomes scarcer. Approaching the limit 

^ a very large lattice is required in order for species 
3 to be viable. In this limit a new, interesting spatial 
organization is observed. Species 3 propagates through 
the lattice in thin and broken wave fronts in constant 
flight from species 2. In the rest of the system the slowly 
growing species 1 and its prey, species 2, is tangled in a 
complex configuration with an enormous mutual perime- 
ter. This spatial organization forms an ever-changing 
labyrinth of narrow pathways in which species 3 propa- 
gates (see Fig. [l]i-f). The more narrow and twisted the 
labyrinth becomes, the longer it will take for species 3 to 
return to a particular location, which gives species 1 more 
time to grow, forming broader pathways. This complex, 
spatial feedback loop stabilizes the configuration. 

In order to describe this spatial self-organization math- 
ematically, we study the probabilities pi , p2 , and ps of a 
random node to be occupied by species 1, 2, or 3, respec- 
tively. Furthermore, we are interested in the perimeter 
Pij between species i and j. That is, the probability of a 
random node and a random of its neighbors to be occu- 
pied by species i and j, respectively. 

Given these perimeters the time evolution of species 
abundances is given by [30] 

Pi = V2P12 - vsPsi^ (1) 

where the equations for p2 and ps follow by cyclic per- 
mutation of the indices 1, 2, and 3. This symmetry also 
holds for all subsequent equations of this article. 
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Figure 1: (Color online) Spatial self-organization in the 
rock-paper-scissors game, a) The three species interact 
cyclically. Species i invades its prey at rate Vi. b-d) Snap- 
shots of the steady state spatial organization of the three 
species when b) All species grow at same rate and L = 300. 

c) Species 3 grows 5 times faster than 1 and 2 and L = 300. 

d) Species 1 grows 1250 times slower than 2 and 3 and 
L — 12800. e-f) zooms of the system in panel d. 



To explain this behavior, one can adopt a mean field 
approximation, where all nodes are linked and spatial 
structure does not exist. Then, the perimeter between 
two species is simply given by the product of species 
abundances pi2 = where tilde (~) denotes that 

the mean field approximation has been applied. If this is 
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Figure 2: (Color online) Species abundancies and 
perimeters for small vi, a) When the growth rate of 
species 1 is decreased, species 3 becomes less abundant. The 
mean field theory correctly predicts the abundancies. b) The 
perimeters between the species are much lower than predicted 
by the mean field theory due to clustering, c) The ratio be- 
tween the predicted mean field perimeters and the observed 
perimeters are equal for all species. This ratio defines x- d) 
When V2, ^ t'i,t'2, the ratio x diverges corresponding to the 
large clustering in Fig. [l]:. When one species grows much 
slower than the others x approaches 5, which gives rise to the 
labyrinthine clustering in Fig. [TJl. 



inserted into ([T]) and the time derivatives are set to zero, 
one obtains the steady state solution 
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In Fig. |2^-b the steady state abundances and perime- 
ters are shown at constant V2 = vs = I and varying 
vi < 1. It is seen that a slow growth rate of species 1 leads 
to a decline in the abundance of species 3 as expected. 
The mean field approximation correctly predicts how the 
abundances of the three species depend on the growth 
rates. However, the mean field approach can not capture 
the spatial organization of the species, and thus it pre- 
dicts perimeters far longer than what is observed in sim- 
ulations (see Fig. The fact that the abundances are 
correctly predicted indicates that the mean field perime- 
ters are proportional to the true, spatial perimeters. In- 
deed, if Q is set to zero for both the spatial and mean 
field system one can derive the relations. 
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Here, we have introduced the ratio x^ defined by how 
much the perimeter between two species is longer in the 
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steady state of the mean field approximation compared 
to the spatial system (see Fig. This new statistical 
measure describes the spatial and dynamical organization 
of the rock-paper-scissors game for varying growth rates. 
The intuition behind x is the following: 

The average time before a node of species 2 is invaded 
by species 1 is given by Ti = ^^^^^ for the spatial system 

and Ti = ^~ in mean field. Therefore, y ~ ^ provides 

1 Vipi2 ' ^ Ti ^ 

a measure for how much longer each species on average 
lives on a node before being invaded, compared to the 
mean field system, i. e. how much the spatial organization 
slows down the dynamics. Furthermore, when x is large 
the perimeters of the spatial system is much smaller than 
in the mean field system, according to Q, so the species 
must have a high degree of clustering. Hence, x gives a 
measure for the clustering of the spatial system. These 
two interpretations are, of course, tightly connected. If 
the average cluster diameters are doubled, each node will 
live for twice as long before being invaded, corresponding 
to increasing x by a factor two. 

How does X depend on the growth rates of the three 
species? In Fig. [2]i this dependency is shown as a func- 
tion of the relative growth rates — and — , with chosen 
to be the fastest growing species. When all growth rates 
are equal we have x ^ 2.5, corresponding to the mod- 
erate amount of clustering observed in Fig. [TJd. When 
species 3 grows much faster than the two other species, 
such that both growth ratios go to zero, x becomes very 
large. This agrees well with the large amount of cluster- 
ing observed in Fig. [T]^. 

When only vi ^ we see from Fig. [T]i and [2]i that 
X approaches a finite value close to 5. In this limit, we 
expect from ([2| that ps ^ while Pi ^ P2 ^ In this 
case, it is limited how much species 3 can cluster. The 
observation of x suggests that clustering only reduces the 
perimeter between species 2 and 3 by a factor 5 compared 
to the mean field system. This sets an upper bound for 
how much species 1 and 2 can cluster. The mean field 
approach predicts a perimeter pi2 = ^, so with x ^ ^ 
equation Q dictates the perimeter in the spatial system 
to be pi2 ^ ^. This agrees well with the 12800 x 12800 
system in Fig. [T]i, where vi = 0.0008, V2 = vs = 1^ and 
Pi2 ^ 0.05, which is also evident from Fig. [SJd. 

While the value of x quantifies the average amount of 
clustering, it does not provide information on the clus- 
ter size distribution. In the case where all species grow 
with the same rate. Fig. [T]3 suggests that clusters have a 
characteristic size. Indeed, Fig. [3] shows that the cluster 
size distribution in this case sharply decreases for clus- 
ters larger than 1000 nodes. When the growth rate of 
species 1 goes to zero, however, species 3 continues to be 
organized in small clusters, but large clusters of species 1 
and 2 become much more likely. The cluster size distri- 
butions of these become exceedingly broad culminating 
in a heavy tail distribution with a cut-off that is set by 
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Figure 3: (Color online) Cluster size distributions, a 

When all species grow at the same rate, all clusters consist of 
less than 5000 nodes. Here vi = V2 = vs = I. b When the 
growth rate of species 1 is decreased, species 3 becomes less 
abundant and large clusters of species 1 and 2 become more 
likely. Here vi ^ 0.5, and V2 = V3 = I. c In the limit ^ 
the cluster size distributions of species 1 and 2 become heavy 
tailed with a cut-off set by the system size. Here vi < 0.007, 
V2=V3 = l. For all plots L 2048. 

the system size. 

An alternative approach that has been applied to de- 
scribe the spatial organisation of the rock-paper-scissors 
game is the pair approximation (STJ [34] . Here the time 
evolution of the perimeters pij are expressed through the 
probabilities p*^ of a random node belonging to species i 
and a random pair of its neighbors belinging to species j 
and /c, respectively. When all growth rates are equal, the 
pair approximation gives the steady state probabilities 

m 

PI = 3, PU = ^, P12 = ^. (5) 

Since (|3| pi2 = ^, the pair approximation predicts x = 
1.5, which is far from the observed value of x = 2.5. This 
further illustrates that the incapability of the pair ap- 
proximation to describe the behavior of the rock-paper- 
scissors game. 
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Discussion — Our results quantitatively describe how 
spatial clustering slows down the dynamics of the rock- 
paper-scissors game, and why this leads to a labyrinthine 
spatial organization in the limit where one species grows 
slowly compared to the two others. An organization that 
includes a new type of excitable fronts that propagate on 
self-organized labyrinthine clusters distributed over many 
length scales. 

In this limit of one slow species, the largest clusters of 
both the slow species and its prey cover a large fraction 
of the system, as seen in Fig. 2d. This consequence of the 
labyrinthine configuration would not be possible in site 
percolation, where each of the large species would need 
to occupy close to 60% of the nodes to percolate [35 . 

Interestingly, the extreme version of the rock-paper- 
scissors ecology with one slow species bear resemblance to 
the forest fire model in a fire-tree- ashes analogy [36H381 . 
The slow species would then be forest, which is burned 
by fire, which is replaced by ashes, from which trees can 
again slowly grow. The main differences from existing 
forest fire models are that in the present system trees 
can only grow in the neighborhood of other trees and fire 
can only be extinguished in the neighborhood of ashes. 

The method of quantifying how much clustering slows 
down the dynamics of a spatial system, compared to the 
mean field approximation, is quite general, and we ex- 
pect it to be applicable on a broad range of dynamical 
systems. In particular, it may be useful predicting the 
spatial organization in predator-prey models, which con- 
tinues to attract much attention within the field of com- 
plex systems [39l [40] . 
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